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ABSTRACT 

We present detailed fits of the spectral energy distributions (SEDs) of four submillimeter (submm) 
galaxies selected by the presence of a gamma-ray burst (GRB) event (GRBs 980703, 000210, 000418 
and 010222). These faint ~ 3 mJy submm emitters at redshift ~ 1 are characterized by an unusual 
combination of long- and short-wavelength properties, namely enhanced submm and/or radio emission 
combined with optical faintness and blue colors. We exclude an active galactic nucleus as the source 
of long-wavelength emission. From the SED fits we conclude that the four galaxies are young (ages 
< 2 Gyr), highly starforming (star formation rates ^ 150 Mq yr^^), low-mass (stellar masses 
~ 10^° Mq) and dusty (dust masses 3 x 10® Mq). Their high dust temperatures (T^ > 45 K) 
indicate that GRB host galaxies are hotter, younger, and less massive counterparts to submm-selected 
galaxies detected so far. Future facilities like Herschel, JCMT/SCUBA-2 and ALMA will test this 
hypothesis enabling measurement of dust temperatures of fainter GRB-selected galaxies. 
Subject headings: galaxies: evolution — galaxies: high-redshift — galaxies: starburst — galaxies: ISM 
— dust, extinction — gamma-ray: burst 



1. INTRODUCTION 

It has been c laimed that submi llimeter (submm) galax- 
ies (SMGs, see IB lain et al.l [200^ for a review) are signif- 
icant contributors to the star fo rmation history of the 
Universe at redshifts z ~ 2 — 3 ([Chapman et al.ll2005f ) 
and have built up a su bstantial fraction of the present- 
day stellar population (|Lillv et al.lll999t ). SMGs are lu- 
minous (with infrared luminosities Lir, ~ 10^^"^'' L©) 
and cold (with mean dust temperature = 36 ± 7 K, 
IChapman et al.ll2005| ). Galaxies with similar luminosi- 
ties but with higher dust temperatures (T^ > 45 K) are 
difficult to detect in the submm with current technology 
due to the fact that the peak of the infrared dust emis- 
sion is shifted out of the sensitive 850 iim band towards 
shorter wavelengths f Blain et al.l 12004 IChapman et al.l 
[2004, 2005) in such galaxies. 

At the other end of the galaxy luminosity function, 
the host galaxies of long-duration gamma-ray bursts 
(GRBs, originating in the collapses o f very massive stars 
at the end of their evolution, e.g. IStanek et al.1 l2003l : 
iHiorth et al.ll200^b[ ) do not have much in common with 
SMGs except for the fact that this type of galaxies is 
also thought to contribute significantl y to, or at least 
trace , the g lobal star formation (e.g. iJakobsson et al.l 
[2005ll2006bh . In contrast to SMGs, GRB hosts are found 
to be blue, subl uminous in the optical (iLe Floc'h et al 



20031: IFruchter et al.ll200 6l and metal-poor ([Fvnbo et al 



20031 ). The majority of them have not been detected 



at mid-infrared (MIR), submm or radio wavelengths 
([Hanlon et al.' '2000^, 'Le Flo c'h et~aLl 120061: iTanvir et all 
[2004; Bcrgcr ct al. 2003; P riddev et al.ll2006f l indicating 
that, as a class, they are not heavily obscured or vio- 
lently star forming galaxies. A low internal dust con- 
tent is consistent with lo w extinction found in the analy- 
sis of optical af terglows ([Stratta et al.ll200^ IChen et al.l 
l2006t iKann et ' al. 2006.) and optical s pectral energy dis- 
tributions (SEDs) of the host galaxies ([Christensen et aT] 
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I2OOI . 

However, four GRB hosts (980703, 000210, 000418 
and 010222) have been firmly detected in submm 
and/ or radio ([Tanvir et all l2'004t iBerger et all l2001bl . 
[2003h providing a somewhat complex picture: as- 
suming that this emission is powered by starbursts, 
the derived star formation rates (SFRs) are of the 
order of a few hundred solar masses per year and 
the amount of dust in these galaxies is significant. 
On the other hand, they exhibit blue colors, low 
extincti on and low ex ti nction - corrected opti c al/UV 
SFRs (iDiorgovski elall [19981: iHolland et al.l mM 
Sokolgv et al. 2001 : ICharv et al l 120021: iGorosabel et al l 



2003allbl: Galama et all 120031 iBerger et all 120031: 



Savaglio et al..,2003i: IChristensen et al .1 120041 : IChen eTall 



20061 iKann et al.l l2006l ). like the majority of known 



GRB hosts. This p uzzling situation was highlighted by 
IBerger et all (|2003l their Figure 3) — these GRB hosts 
are much fainter in the optical than the prototypical 
Ultra Luminous Infrared Galaxy (ULIRG) Arp 220 
and have much bluer spectra, but are more luminous 
in su bmm and radio. I n fact, no t emplate SED model 
(e^ lSilva et~aI1 [l998l : iDale et al.l [20011 : iDale fc Heloiil 



2002f) could giv e consistenc y with the data (see also 
Michalowski fc Hiorth 2007/lMichalowskil [20061 ) . More- 



over, none of th e se ho st s have been detecte d at 2 4 /zm 
(|Le Floc'h et al.l 120061: ICastro Ceron et al.l l2006l . Le 

Floc'h, private communication). 

The location of the GRB events within their hosts 
ad ds further cornplexity to this picture. As found 
by IFruchter et al.l (|2006D . GRBs trace the location of 
the brightest res t-UV parts of their hosts (see also 
iBloom et al.ll2002[ ). If the majority of the star formation 
in the hosts of GRBs 980703, 000210, 000418 and 010222 
was hidden by dust, then they should preferentially be 
found in obscured (UV-dim) parts of their hosts (as long 
as GRBs trace star formation), contrary to observations. 

In this paper we investigate this seeming discrepancy 
between short- and long-wavelength data through stel- 
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lar population model SED fitting. In particular we dis- 
cuss the possibility that these submm-bright GRB hosts 
may represent the long-sought hotter (and less luminous) 
counterparts of SMGs. In Section [2] we describe the 
model and the fitting procedure, in Section [3] we show 
the results, discuss their implication in Section [4l and 
finally in Section [5] conclude this work. We use a cosmo- 
logical model with Ho ~ 70 km s~^ Mpc"""^, J^a — 0.7 
and i^M — 0.3. 



2. GRASIL SED MODELING 

In order to model the SEDs of GRB hosts we us ed the 
GRASIL 1 software developed bv lSilva et all (|1998D . It is 
a numerical code that calculates the spectrum of a galaxy 
by means of a two-dimensional radiative transfer method, 
applied to photons produced by a stellar population, and 
reprocessed by dust. This model is self-consistent in that 
it fulfills the principle of energy conservation between the 
energy absorbed by dust in the UV/optical wavelengths 
and the energy re-emitted in the infrared. Two extinc- 
tion media are implemented, dense star forming molec- 
ular clouds (MCs, applied only to the youngest stellar 
population) and diffuse cirrus. The dust is composed of 
small grains (not in thermal equilibrium with radiation 
and hence fluctuating in temperature), big grains (sili- 
cates and graphites) and polycyclic aromatic hydrocar- 
bon (PAH) molecules. The emission of grains with given 
size is assumed to be a grey-body, so the composite spec- 
trum is a sum of all grey-bodies with different tempera- 
tures. A galaxy is assumed to be an axially symmetric 
system. Different geometries of the stellar and diffuse 
dust distributions are allowed but were not used here. 
The SFR is assumed to be proportional to the available 
gas content following the Schmidt-law. On top of this 
smooth SFR history a violent starburst epoch is added. 
Star formation is unevenly distributed throughout the 
galaxy in MCs. The SFR given as a GRASIL output is 
the sum of the SFRs of each MC. 

We gathered photometr ic data f r om the litera- 
ture at optical dSok olov et al.| 12 0011: IVreeswiik et al.l 
Il999t iGorosabel et al.ll2 003a b: Gal ama et al.ll2003D. in - 
frared (ICastro Ceron et al. 2006; Le Floc'h et al.ll2006D . 
submm (iTanvir et a l. 2004; Berger et al. 2003 ) and radio 



(iBereer et al.ll2001bl I2003t iFrail et al.ii2002i : iSagar et all 
|200l|) wavelengths. For the host of GRB 000210 we 
performe d the photometry on th e arch ival Spitzer im- 
ages (see lCastro Ceron et al.ll2006ll200 7l for a decription 
of the procedure) and obtained the following fluxes: at 
4.5 /im: 3.3 ± 2.1 /iJy; at 8.0 /xm: 15.0 ± 5.1 /iJy and at 
24 ^m: < 31.5 ^Jy (3cr). 

We performed SED modeling investigating a wide 
range of the following GRASIL parameters: age of the 
galaxy (defined as the time since the beginning of its evo- 
lution when the stellar population started to build up), 
dust-to-gas ratio, and mass of gas converted into stars 
during the current starburst episode (lasting 50 Myr). 
We used a clearing time for MCs, tesc — 50 Myr (see 
iPanuzzo et al.|[2007l . for discussion of this parameter). 



The best fits^ are shown in Figure [T] and the results 
for each parameter listed in Table [T] together with sev- 
eral properties of the galaxies derived from the SEDs. 
SFRs, stellar masses, dust masses and extinction in the 
MCs are given as output from GRASIL. Infrared lumi- 
nosities were obtained by integrating the SEDs over the 
range of 8 — 1000 fim. Dust temperatures were estimated 
by fitting a grey-body curve to t he part of the SEDs 
near the dust peak (^ 100 fim, e.g. lYun k Cariilil [2b021. 
Finally, the average extinction outside MCs was calcu- 
lated as: Av — 2.5 log(V^-band starlight exti nguished by 
MCs only / F-band starhght observed); see iSilva et al.l 
(|1998f) . This parameter describes the extinction averaged 
throughout the galaxy as opposed to the line-of-sight ex- 
tinction derived from optical GRB afterglows. 

We checked the consistency of our resu l ts with those 
reported in the literature (Berger et al.l | 2001aL 1 2003 : 
Biornsson e t all 120021; ICastro Ceron et al.' '2000, '2007; 



3. RESULTS 

^ |http: //web. pd. astro . It/granato/grasil/grasll .html | 



Chen et al . 2006; Christensen et al. 2004; Galama ct al] 
2003; Gorosabel et al. 2003a,b; Le Floc'h et al. 200|; 
Michal owski 2006; Michalows ki fc Hiorth ,20071: 
Sokolov et all l200ll: IStratta et all l2004l : iTakagi et al l 
2004 see Appendix [A|. We generally found good 
agreement except for the following points. Our esti- 
mate of the age for GRB 980703 (which has the most 
significant old stellar c omponent) is cons i derabl y larger 
than that derived by IChristensen et all ()2004f ). The 
discrepancy arises because we calculated the time from 
the beginning of the galaxy evolution, not the beginning 
of the starburst. Our SFRs for GRB 980703 and 010222 
are inconsisten t with the reported upper limits of 
iLe Floc'h et al.l (|2006f) . However our analysis, based on 
the full SED rather than only 24 /xm d atapoints, seems 
to be more reliable as admitted by iLe Floc'h et al.l 
(|2006h who assigned a factor of ^, 5 error to their 
estimates (see also the results of ICastro Ceron et al.l 
120061 based on the same data). iBerger et all (120031 ) 
obtained systematically higher SFRs by a factor of 
2-5 based on submm alone. This may be the effect of 
uncertainty of the dust pro p erties (we co uld reproduce 
the re sults of IBerger et al.l (j2003D using lYun fc Carillil 
(j2002|) . but when we used a temperature 10 K higher 
and a n emissivity index 0.35 lower than IBerger et al.l 
(I2OOI did, we obtained values consistent with our 
results reported here) as well as of the contribution of 
older stellar populations to the submm fluxes ( see for 
example the discussion in IVlahakis et"al] l2007l ). that 
leads to an overestimation of SFRs using only submm 
data. For the hosts of GRB 000210, 000418 and 010222, 
ITakagi et"an (|2004f ) predicted much higher SFRs, of the 
order of a few thousand M0 yr^^ and stellar masses 
of ~10"Mq. Such a high SFR is unlikely to be 
necessary to explain the submm emission of GRB hosts 
and hence the accumulated stellar masses seem to be 
too high as well. We derived much higher extinction 
for th e host of GRB 000210 than IChristensen et all 
dlOOl), but our model predicts that the majority of 
the extinction in the V^-band has a grey nature, which 
would be undetectable in any reddening measurements 
(we found E{B — V) = 0.07, which is consistent with 
their results if one assumes the Galactic extinction curve 

^ The SED fits can be downloaded from 
|http : / /archive . dark- cosmology . dk | 
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Fig. 1. — SEDs of GRB hosts. Solid lines: the yo ung starburst galaxy models calculated using GRASIL and consistent with the data. 
Dotted lines: Arp 220 model (from lSilva et allllQQSI) . Squares: detections with errors, in most cases, smaller than the size of the symbols. 
Arrows: 3cr upper limits (values marked at the base). Hashed columns mark the wavelengths corresponding to optical/near-infrared filters, 
MIR Spitzer filters, SCUBA submm bands and the radio domain. It is apparent that the ULIRG Arp 220 model underpredicts the radio 
and submm fluxes of GRB hosts while fitted to the optical data. The optical color of this model is also too red compared to GRB hosts. 
In contrast, the emission of young starbursts is blue in optical and strong in submm/radio at the same time. In the case where we lack 
Spitzer data the MIR part of the SED is unconstrained due to the possible absorption or emission of PAH and silicate features. 



slope Rv = 3.1). Moreover, it is possible that a GRB 
event destroys the dust along the line of sight, so study 
of an afterglow results in low dust content. 

The determination of the infrared luminosity suffers 
from systematic uncertainties depending on the choice of 
the SED template. Our approach of using all the op- 
tical, submm and radio data to constrain the shape of 
the SED results in a moderate systematic e rror i n lu- 
minosity. Using the templates of iDale et al.l ()2001[ ) and 
I Dale fc Heloul ()2002[ ) fitted to long-wavelength data we 
obtained values only 30% l ower. Sim i lar an alysis on a 
bigger sample of galaxies led lBell et al.l ()2007D to the con- 
clusion that the systematic uncertainty of infrared lumi- 
nosity is usually less than a factor of ~ 2. The arbitrary 
choice of the^alpctcr (1955) IMF with the cutoffs of 0.15 
and 120 Mq introduces a systematic error of a factor of 
^ 2 in the deter m ination of t he ste llar masses and SFRs 
(|Erb et al.ll2006l ). iBell et al.l ()2007f ) have also found that 
random errors in stellar mass are less than a factor of 



~ 2. Dust mass estima tes are uncertain up to a factor of 
a few (jSilva et al.lll998D . Estimates of dust temperatures 
based o n submm and radio a lone have uncertainties of 
- 10 K (jChapman et al.l[2003h . 

4. DISCUSSION 

4.1. Solving the puzzle 

The key property of GRB hosts that explains their 
blue colors and enhanced submm/radio emission is their 
young ages (Table [T]). On one hand the majority of the 
stars still reside in dense MCs, so a significant part of the 
energy is absorbed and re-emitted. This increases the 
dust emission. On the other hand there are lots of young, 
hot, blue stars in such a galaxy, because they have not 
finished their lives yet. Hence the total optical spectrum 
of the galaxy is blue. GRBs may indeed reside in or close 
to MCs (possibly causing hydrogen ionization and dust 
sublimation along the line of sight. IWatson et "all l2007l ) 
and it was found that gas column densities derived from 
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TABLE 1 

GRASIL PARAMETERS YIELDING SEDS CONSISTENT WITH THE DATA AND CHARACTERISTICS OF THE GALAXIES DERIVED FROM THE SED 

MODELING. 



GRB 


z 


Age 


dust / gas 






SFR 


M* 




Ti 


/(MO 


A- 


host 




(Gyr) 


(10-2) 


(109 Mo) 


(1012Lq) 


(M© yr-i) 


(109 Mq) 


(109 Mq) 


(K) 




(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


(8) 


(9) 


(10) 


(11) 


(12) 


980703 


0.97 


2 


0.3 


2.3 


1.9 


212 


21.0 


0.15 


44 


54 


0.42 


000210 


0.85 


0.19 


2.4 


1.2 


0.9 


179 


10.1 


0.33 


45 


27 


1.01 


000418 


1.12 


0.14 


1 


5.0 


4.6 


288 


18.7 


0.82 


50 


31 


0.70 


010222 


1.48 


0.09 


0.7 


5.9 


4.3 


278 


13.5 


0.32 


51 


25 


0.58 


Arp 220 


0.018 


13 


1 


25 


1.16 


580 


153.2 


0.30 


55 


63 


0.59 


NGC 6946 


0.00016 


13 


1 




0.05 


6 


104.1 


0.10 


29 


36 


0.06 



Not e. — (1) GRB name, (2) redshifts from lDiorgovski et al.l l lT99g) . |PFo et al.l | |2003) , IBloom et all l l200l) , IJha et al.l 120011) , ISoifer et all 
111 9841) and lHuchra et aTl II1999I) . (3) age of the galaxy defined as the time since the beginning of its evolution, (4) dust-to-gas ratio, (5) mass 
of gas converted into stars during starburst, (6) total 8 — 1000 /jm infrared luminosity (7) total star formation rate for 0.15 — 120 Mq stars 
averaged over the last 50 Myr, (8) stellar mass, (9) dust mass, (10) dust temperature, (11) extinction in molecular cloud (MC) at 1 /im 
measured from its center, (12) average extinction of stars outside MCs at 0.55 /xm. Values corresponding to Arp 220 (ULIRG) and NGC 
6946 (spiral) are given for comparison. 



X-ray afterglows in a sample of 8 GRBs (including GRB 
980703, discussed here) were in the range corresponding 
to the column densities of giant MCs in the Milky Way 
(jGalama fc Wiieril2001h . A similar conclusion for high- 
redshift GRBs was recently drawn by Jakobsso n et al.l 
(|2006ari W means of modeling Lyg absorption features 
(see also [Castro -Tirado et al."1999l: 'Hiorth et al."l2003j: 
Savaglio et"aL|[2 003: Stratta ct al. 2004: Vrccswiik ct aj] 
2004 IChen et al. 2005, 2006; Wats on et all 12001 
Campana et al., , 2007; Prochaska ct 'aLl l2007allbl : 
Ruiz-Velasco et al 1 120071 for a discussion of similar 



results). 

Since the stars formed during the starburst do not 
dominate the stellar mass of the GRB hosts discussed 
here (only from 10 to 40% of these masses have been 
formed during the on-going starburst epoch — com- 
pare columns 5 and 8 of Table [1]) and are still embed- 
ded in MCs providing st rong extinction (so-c alled age- 
dependent extinction, see'Panuz zo et ani2007f l. it is ap- 
parent that the optical/UV light is dominated by some- 
what older stars, which are already outside MCs and 
suf fer only moderate extinc tion (column 12 in Table [TJ 
see iPlante fc Sauvagj |2002| . for an example of a totally 
obscured, metal-poor star-forming region). This is why 
optical measurements of extinction resulted in low values 
suggesting low dust content whereas enhanced submm 
emission is consistent with being emitted by dusty galax- 
ie s. In light of this, we prov ide support for the hypothesis 
of iGorosabel et aTl ()2003 a') claiming that the optical and 
submm emission of the GRB 000210 host are dominated 
by different popul ations of stars (the s ame explanation 
was proposed by iGoldader et al.l I2002L for a sample of 
dusty ULIRGs with UV colors bluer than expected for 
starburst galaxies). 

The l ocation of GRBs in th e brightest UV parts of the 
galaxy (jPruchter et al.l l2006') is also consistent with our 
model. Although the GRBs discussed here trace regions 
of obscured rather than non-obscured star formation (be- 
cause the majority of star formation is obscured and 
under the assumption that GRBs trace star formation), 
these regions are not spatially distinct in the galaxy. Re- 
gions of obscured star formation evolve into non-obscured 
regions by destroying the MCs without changing their lo- 
cation. Unless individual MCs can be resolved, obscured 



star formation within them and less-obscured star for- 
mation on their outskirts cannot be spatially separated. 

Our results are based on the assumption that the de- 
tected submm/radio sources are indeed related to GRB 
hosts. It is however possible that the emission comes 
from unrelated sources falli ng into the c o arse beam o f 
SCUBA (15") a s noted by iSrnith et all (l20?)ll . [20051), 
IGorosabel efall and iLe Floc'h et al.l ([2006). 

The most suspicious case is the GRB 010222 host which 
is undetected at 24 /im, but accompanied by several spa- 
tially close MIR-b right galaxies that cou ld dominate the 
submm emission (|Le Floc'h et al.ll2006[ ). On the other 
hand, the nondetections of GRB hosts in the MIR are 
easily explained by silicate absorption features and the 
steep infrared spectrum of the galaxies (see Figure [1]). 
The situation is less severe for GRB 980703 and GRB 
000418, both detected by VLA in the radio, for which 
accurate astrometry decreases the chance of confusion. 

4.2. The nature of the GRB hosts 

From the results presented in Table [T] a common 
characteristic of submm GRB hosts emerges (although 
our sample is too small to draw a very robust conclu- 
sion). For all four galaxies we obtained young ages 
2 Gyr) and relatively low stellar masses (< 2 X 
10^°M f7i). These are usual properties among GRB 
hosts (IChristensen et al l |2004 ISollerman et afl 120051 : 
ICastro Ceron et al.ll2006l : ISavagho et al.ll2006l ) .Galaxies 
with such stellar mass are believed to domi nate the star 
formation of the Universe at redshifts z <1 (jZheng et al.l 
2007; Bell ct al. 2007!). 

What is special about the four galaxies discussed here 
is that they are highly star formi ng (as oppo s ed to 
the other known GRB hosts, see iBerger et al.l 120031: 
ILe Floc'h et~all 120061 : iTanvir et all 120041 ) andlmve high 
dust content. Hence they can be representative only 
of the bright end of the infrared luminosity function of 
GRB hosts. More sensitive observations in submm (by 
JCMT/SCUBA-2 and ALMA) and far-infrared (FIR, by 
Herschel) are therefore necessary to detect fainter hosts 
and test this hypothesis. 

The optic al afterglow of GRB 000418 was excep- 
tionally red (|Klose et al.l [2000) while no optical after- 
glow was detected for GRB 000210 in spite of deep 
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searches (jPiro et alj |2002[ ) (its X-ray luminosity was 
high enough to place it at the borderline of the cat- 
of dark GRBs, iJakobsson eFan 12004 iRol et all 
both hinting at significant obscuration in the 
hosts. However the "darkness" of GRBs cannot eas- 
ily be linked with obscuration in the hosts, because 
the remaining three members of our submm/radio 
sample were associate d with br ight optical afte rglows 
(iBerger et al.ll2 001a: Bl oom et a l. 1998; Biornsso n et al.l 



rg( 

200a iCastro-T irado et al.' 1999; Galama e t all 12005 
Holland et aLri2G0I; Klose et al. 2000; Vreeswiik et all 



19991). Moreover, iBarnard et all (|2003l ) did not detect 



submm emission from the hosts of three GRBs classified 
as dark. This issue should be addressed by targeting a 
more significant sample of dark GRBs in the submm and 
radio. 

4.3. Rejection of an AGN contribution 

In general, GRB hosts are typically starburst galaxies 
and this selection makes the presence of an AGN compo- 
nent quite unlikely. Here we present additional indica- 
tions that the long-wavelength emission from GRB hosts 
discussed here is not dominated by AGNs. 

In order to assess the probability that there is an AGN 
component in our sample of GRB hosts we compared 
th eir FIR and r adio l uminosities using the q coefficient 
of iHelou et all (fl98l : q = log(FIR/3.75 • W^/Sia)- 
The FIR luminosity was integrated in the range 42.5- 
122.5 //m. The rest 1.4 GHz luminosity {Si 4) was calcu- 
lated from the observed 1.4 GHz (or 4.86 GHz for GRB 
010222) flux assuming a steep radio slope a = —0.75 to 
obtain a robust lower limit on q. (Assuming a shallower 
slope, the expected rest 1.4 GHz radio luminosity would 
be even lower). The resulting q values are 2.28, > 1.93, 
2.41 and > 2.06 for GRBs 980703, 000210, 000418 and 
010222, respectively. It is known that starburst galaxies 
follow the so-ca lled FIR-radio co rrelation with a mean 
Q = 2 .21 ± 0.14 (|Helou et al.lll985n or 2.3 ± 0.2 (|Condonl 
Il992f ). Hence the GRB hosts are consistent with be- 
ing starburst-dominated galaxies. Emi ssion dominated 
by radio-loud AGN would show q < 2 (|Yun e t al. 2001; 
IChanial et aI.H2007l : lYang et al.ll2007t ) which would only 
be the case for GRB 000210 if the actual radio flux is 
just below the 3a upper limit. We therefore conclude 
that the FIR-radio correlation shows that the emission 
from GRB hosts is probably not dominated by radio-loud 
AGNs. Even if a non-dominating AGN is present its con- 
tribution to the radio and submm emission is insignifl- 
cant since the hosts fulfill the FIR-radio correlation. It 
would be rather unlikely that AGN-domianted emission 
is coincidentally consistent with tight FIR-radio correla- 
tion. 

This is supporte d by several other diagnostics. 
iBerger et al.l ()2001bD detected no variability in the ra- 
dio flux of GRB 980703 host over 650 days contrary to 
what would be expected if an AGN dominates its ra- 
dio emission. Moreover, GRB hosts have optical spectra 
typical for starbursts, not AGN s (i.e. no high-ion i zation 
AGN lines have been found. iDiorgoyski et al.l [19981: 
Berger et ai]|2001bt [Bloom et al.l 120031: iProchaska et all 



200- 



This excludes non-obscured AGNs in our sam- 
ple. Moreover, GRB hosts are dwarf galaxies and the 
fraction of AGNs (without optical high ionization lines) 
in a galaxy sample with stellar mass of M^, ~ 10^" Mq 



is negli gible (< 0.01% fr om Figure 2 of [Best et all 120051 : 
see also lWoo et all[2005l) . AGNs are also more luminous 
than GRB hosts as is shown in Figure [H In summary, 
AGNs cannot dominate the emission of the GRB hosts 
discussed here unless they are atypically small, obscured 
and radio-quiet. Since submm- and radio- faint hosts are 
even less probably connected with AGN activity (because 
it would require even smaller AGNs) we conclude that 
GRB hosts in general are unlikely powered by AGNs. 

4.4. The general picture of dust properties 

In Figure [2| we compare the total infrared luminosities 
and dust temperatures of GRB hosts with well-studied 
galaxies both local and at high-z. We included the four 
hosts studied here as well as the host of GRB 030115 
which has an upper limit t o the tempera t ure of 50 K. 
This limit was derived by iPriddev et al.l (|2006[ ) using 
SED modeling of optical and near-infrared data, which 
led to the estimation of the infrared luminosity via the 
UV-slope method. This, in turn, allowed them to ex- 
clude a low dust temperature, because cold dust would 
be inconsistent with the non-detection in the submm. 

GRB hosts seem to overlap with intermediate-z 
ULIRGs from lYang et all ([2007). This is not very sur- 
prising since both galaxy classes have ULIRG charac- 
teristics. GRB hosts are however much more distant. 
As opposed to a mean redshift of 0.37 for intermediate-z 
ULIRGs, our sample has a r nean redshift of 1.1 comp ared 
to 2.8 for GRBs in general (IJakobsson et al.ll2006bl ). 2.2 



for SMGs (jChapman e t all 120051) and 2.1 for opt ically 
faint radio galaxies (OFRGs, iChapmaii "eral1[200l . 

The GRB hosts discussed here occupy th e same region 
of Fig ure [2| as the brightest starbursts from lTavlor et all 
(|2005f) . Moreover, dust masses of GRB hosts (Tabled]) 
are close to the upper boundary for starburst galaxies 
derived by these authors. 

As noted above, the GRB hosts discussed here must 
form the bright end of the infrared luminosity function 
of GRB hosts. The remaining members of the sam- 
ple are still undetected at long wavelengths making it 
impossible to measure luminosities and dust tempera- 
tures. One can speculate from Figure [D that they may 
follow the "starburst" sequence having similar tempera- 
tures (T ^ 40 — 50 K) and infrared luminosities in the 
range 10^ — 10^^ Lq. Herschel should be able to detect 
these sources in the FIR. There is also the possibility 
that we can detect even brighter (but rare) G RB hosts — 
the co unterparts of the brightest ULIRGs of I Yang et all 
(|2007D with luminosities ^ 10^^ Lq and temperatures 
~ 60 K. These should be bright enough to be searched 
for by sensitive submm instruments (JCMT/SCUBA-2, 
ALMA). 

It is known that GRB hosts are much bluer than 
mass ive , star forming SMGs (co mp are I Christensen et al.l 
200I ILe Floc'h et all [2001 IBerger et all 120031 : 
Small et al. I 120041) . From Figure [D it seems that 
they are also hotter than SMGs with the same luminos- 
ity (or dimmer than SMGs with the same temperature). 
This gives a hint that GRB events may pinpoint a 
population of ULIRGs at high redshifts with dust hotter 
than in typical SMGs. The search for such galaxies is 
important because they likely contribute to the star 
formation history at the same level as SMGs. High 
dust temperatures of GRB hosts were hypothesized 
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Fig. 2. — Dust temperature as a function of inf rared (8 — 1000 lim) luminosity. GRB hosts discussed here (fill ed squares) and GRB 
030115 (arro w indicating lower limit on temperature, IPriddeY. et al."200C?) are compared with submm galaxies (pfaaea. IChapma n et al.ll2005l : 
IKneib et al.| [2005. large symbol denotes a hot lenscd galaxy found behind a cluster), optically faint ra dio ga laxies (triangles, note that 
symbols i ndica te lower limits on temperature, Chapman ct al. 2004), intermediate-2 ULIRGs (crosses, [Yang ct al. 2007), local ULIRGs 
(asterisks. Sol omon e t al. 1997), local LIRGs (dots, Dunne ct al. 2000), low-z starburst galaxies (diamonds, Taylor ct al. 2005), low-z spirals 
and cirrus galaxies (stars, S tevens et al. 2005 : Taylor ct al. 2005) and millimeter-selected radio-quite AGNs (circles, Be nford et al.. ,19991 ; 
iBeclcn ct al. 2006). GRB hosts seem to reflect the properties of intermediate-z ULIRGs, the bright end of starburst galaxies and optically 
faint radio galaxies — the candidates for hotter counterparts of SMGs. Moreover, the hot, faint SMG found behind the cluster (large plus) 
falls in the same region as GRB hosts. 

the former galaxies. If they are compact then the steUar 
population suffers strong extinction. This would lead to 
redder colors, but it is hkely (see Section B?T|) that this 
extinction is so strong that very young stars are totally 
obscured, so optical light is dominated by relatively less 
obscured stars outside star forming regions leading to 
blue colors. 

We note that the majority of the galaxies shown in 
Figure[2]also have higher dust temperatures compared to 
SMGs. However, all of them are local galaxies, so cannot 
be considered as counterparts of high-redshift SMGs and 
their submm emission has been detected only because of 
their proximity. 

GRB hosts may be consistent with a population of 
OFRGs having similar infrared luminosities and (likely) 
temperat ures. Altho ugh the majority of OFRGs lie at 
z ~ 2 (C hapman et al. 2004), some of them are within 
the redshift range of the GRB hosts discussed here. 



by IBarnard et all (|2003 h and ' Tanvir et aP 
possible way to explain their faintness in submm. Here 
we provide evidence t hat this is the case. Moreover, 
iTrentham et al.l ()2002l ) suggested that GRB hosts may 
be low-luminosity SMGs. Indeed, very faint sub-mJy 
SMGs magnified by clusters of galaxies are found to be 
hotter than those found in blind surveys (limited by 
confusion to ~ 2 mJy): t he z ^ 2.5 S MG b ehind the 
cluster A 2218 found by iKneib et all ()2005l ) has very 
similar dust properties to the GRB hosts discussed here 
(Trf~50, Lir-IQI^Lq). 

Hotter dust temperatures indicate that the star forma- 
tion sit es in GRB hosts are more compact th an those in 
SMGs (iChanial et al.l[2007l : lYang et al.l[2007l) . It is con- 
sistent with the fact that GRB hosts have higher specific 
SFRs (per unit mass) than SMGs (jCastro Ceron et al.l 
[2P06). The blue optical colors of GRB hosts compared 
to SMGs can also be explained by the compactness of 



GRB-selected submillimeter galaxies 



OFRGs have been suggested to be hotter counterparts 
of SMGs (Ch apman et alJ [20?)1 . so the same may be 
true for GRB hosts. 

Indeed SMG samples are clearly biased against the 
high-Trf-lo w-L galaxies (upper-left part of Figure [21 see 
iChapnian et al. 200a) . The lack of SMGs in low-Trf- 
high-L (lower-right corner) is probably real, because 
such sources would b e detected if present (see Figure 
2 of iBlain et alJ I2004L for a discussion of selection ef- 
fects). Similarly, the lack of very luminous sources with 
L > 10^'' Lq is probably real. If such powerful sources 
exist, they are very rare and do not contribute to the 
presented sample. Taking into account all the galaxies 
shown in Figure [H the lack of high-T^-low-L galaxies 
is probably not a selection effect (at low redshift) be- 
cause the hotter counterparts of 10^° Lq galaxies should 
be easily detected in MIR and FIR. 

Therefore the apparent trend (also seen in our GRB 
host sample) that more luminous galaxies have higher 
dust temperatures is a real effect called t he luminosity- 
temp erature relat ion (ISoifer et al.l Il987t IChanial et al.1 
120071: lYang et all 1200711 The spread has been inter- 
preted as variation in the size of a star forming region 
(|Chanial et al.ll2"007f ) or in the dust emissivity in dex and 
amount of dust in each galaxy (|Yang et al.ll2007t ) . Galax- 
ies with large dust content tend to occupy the lower-right 
corner of Figure [2] whereas those with low dust content 
occupy the upper-left corner. The GRB hosts discussed 
here with Md ~ 10^ Mq are placed near the average of 
all the galaxies shown in Figure [2l The rest of the popu- 
lation is probably aligned along the diagonal of Figure [2] 
fulfilling the luminosity-temperature relation. It is ex- 
pected that some of the IR-faint GRB hosts have much 
lower dust content so are possibly located in the lower- 
L — higher- regime. 
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5. CONCLUSIONS 

The short- and long-wavelength properties of the host 
galaxies of GRBs 980703, 000210, 000418 and 010222 
can be linked assuming that they are very young and 
powerful starbursts. We conclude that they are galax- 
ies with the highest star formation rates among known 
GRB hosts, but their optical properties, starburst na- 
ture, stellar masses and ages are not distinctive. We also 
found that AGNs are probably not responsible for boost- 
ing their long-wavelength emission. 

We have shown that GRB host galaxies at cosmological 
redshifts may constitute a population of hot submillime- 
ter galaxies. This hypothesis makes GRB hosts of special 
interest, placing them in the same category as optically 
faint radio galaxies, and should be confirmed by future 
sensitive long-wavelength observations. Future instru- 
ments (Herschel, JCMT/SCUBA-2 and ALMA) wiU be 
able to build up a statistically significant sample of GRB 
hosts with accurately measured infrared luminosities and 
temperatures. 
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APPENDIX 

COMPARISON OF OUR RESULTS WITH THE LITERATURE 

In the following appendix we compare our results with those found by other authors applying different methods, 
to show that our modeling in most cases gives consistent values, but also provides additional galaxy properties that 
cannot be inferred from previous approaches. We stress that some inconsistency with methods based on optical/UV 
is expected because optical/UV light traces only a minor (i.e. unobscured) portion of the bolometric luminosity of the 
galaxies. Some authors used different IMFs, but the necessary correction to total SFRs and stellar masses is not larger 
than 15%, and does not affect the conclusions of the comparison. 

Ages 

Based on galaxy SEDs fitted to the optical data only,'Gorosabel et al.' (' 2003al lbl) and lChristensen et al.l (|2004j) derived 
the ages of the starbursts in the hosts of GRBs 980703, 000210 and 000418. Our estimates are considerably larger only 
for GRB 980703. The discrepa ncy is because we ca lculated the time from the beginning of the galaxy evolution, not 
the beginning of the starburst. ISokolov et al.1 (|2001h derived ages of both old stellar populations and starbursts. Our 
estimation for GRB 980703 agrees with the age of the old component, which is conceptually closer to our definition 
of the galaxy age. Ages derived bv iTakagi et all ()2004l ) for GRB 000210, 000418 and 010222 agree with our results 
within a factor of a few. 

Star formation rates 

See iMichalowski fc Hjorthl (|2007f l for details of the comparison between SFR estimates. Si nce the starburst is still 
hidden in MCs and its o ptical light is extin guished, SFRs derived from optical indicators (|Christensen et all 12004 
iGorosabel et al.|[2003ai rbl: [Berger et al.ll2003|) are two o rders of magnitude lower than our estimates. Our results are 
consistent with radio-derived SFRs (iBerger et al.ll2003f ). This is because the calibration of SFR to radio flux r equires 
the prior assumption of only two parameters (a normalization factor and a spectral inde x. seelVun fc Car illi 2002*). 
which are relatively well constrained. O ur values also agree w ith the upper limits derived bv lCastro Ceron et al. (20061 ) 
using the template of Arp 220. Finally, iBerger et al.l ()2003[ ) obtained systematically higher SFRs by a factor of 2 — 5 
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based on submm alone. We have checked that the SFRs derived fr om our SEP models using the total infrared emission 
(H enmcutdllQQSl ) are consistent with our values derived here (see iMichalowski fc Hiorthll2007t l. 

Stellar masses 

Our results for the GRB 980703 host agree with th e stell ar mass derived by iCastro Ceron et alj (l2006f) arid with 
both stellar and burst masses derived bv lSokolov et al.l (|2001|) . The stellar mass derived bv lCastro Ceron et al.l ()2007D 
for GRB 000210 is also in agreement with our value. 

Dust properties 

Our value of th e dust mass for GRB 980703 host agrees within a factor of 1.5 with the one derived by 
ICastro Ceron et al.l (gOp). Dust masses derived bv iTakagi etld] (pOP) for GRB 000210, 000418 and 010222 agree 
with our results within a factor o f a few. We have checked that for th ese three hosts, values of dust mass es computed 
direc tly from submm fluxes (using lHildebrandl[l98llTavlor etalMo^ agree with those reported here (see iMichalowskil 
[20061) . 

The near- infrared extinction derived for MCs (colurnn 11 i n Table [H) is within the typical values found in observations 
of compact star forming regions (e.g. 'Scovillc et al.' 1998; 'Plante & bauvage '2002'; 'Hunt et alj I2005D and numerical 
simulations (e.g. [indebctouw et al. 2006; Goicoechea & Lc Bourlot 2007; Panuzzo et al. 2007). Our values of average 
ex tinction outside MCs , Ay (column 1 2 in Table [H), are c onsistent with those derived from optical host SED modeling 
bv lSokolov erall (|2001|. for 980703 ) and lChristensen et al.l (12004 for 980703. 000210, 00418) and fro r a opti cal afterglow 
modehng bv iBerger et al! (l2001al for 000418). iBiornsson et al.] (|2002l for 010222), IStratta et all (|2004 for 980703, 



010222) and lChen et all ()2006l . for 980703), except for the host of GRB 000210 for which we predict higher extinction, 
but with grey nature undetectable in any reddening measurements To the best of our knowledge we present the first 
Ay determination from the host galaxy SED fitting for GRB 010222. 
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